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r~H Recently, a method based on stochastic quantization has been proposed to compute the 

^ Casimir force and its fluctuations in arbitrary geometries. It relies on the spectral de- 

composition of the Laplacian operator in the given geometry. Both quantum and thermal 
^N) fluctuations are considered. Here we use such method to compute the Casimir force on 

the plates of a flnite piston of arbitrary cross section. Asymptotic expressions valid at low 
C " } and high temperatures and short and long distances are obtained. The case of a piston 

with triangular cross section is analyzed in detail. The regularization of the divergent 
stress tensor is described. 
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1. Introduction 



H.G. Casimir ^ predicted in 1948 the existence of an attactive force between two 
perfect parallel conducting plates due to quantum fluctuations of the vacuum 1^. Ad- 
vanced experimental setups allow to measure such forces for geometries other than 
the two parallel plates However, techniques to calculate forces for complicated 
geometries, beyond the usual ones with high degrees of symmetry, were scarce until 
Ref. m introduced a method to calculate electromagnetic (EM) Casimir forces. It is 
based on a multiscattering calculation and has been successfully applied to many 
configurations, such as plates, cylinders, spheres, wedges, etc. ^ 

Casimir forces have their origin in fluctuations of EM fields. Already Lifshitz in 
1956 ^ use the concept of fluctuations (in this case of the electric currents inside 
a metal) to compute such forces. In Ref. 7 a new method to compute Casimir 
forces has been derived by us. It is based on the stochastic quantization approach 
developed by Parisi and Wu ^. A Langevin equation for a given field subjected to 
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thermal fluctuations, is constructed in a non-physical, pseudo-time, such that it 
reproduces the probability distribution for the field in a thermal bath. Using the 
formalism developed in Ref. Elto compute fluctuation- induced forces starting from a 
stochastic equation, the Casimir force of the electromagnetic field can be obtained. 
The resulting method has the advantage of computing the Casimir force directly 
and can be applied to compute torques and stresses over extended bodies as well. 
Also, the force fluctuations are obtained properly regularized and it was shown that 
the force variance on a piston geometry is finite, and twice the square of the force 

In this manuscript we apply this method to the computation of the force on pis- 
tons of arbitrary cross section. The text is organized as follows. Section [2] presents a 
summary of the procedure to obtain the Casimir force using the stochastic quanti- 
zation method. An analysis of the divergence of the stress tensor is performed and 
a force regularization method is proposed. In Section|3]the method is applied to the 
case of the electromagnetic Casimir force on pistons. The decomposition in trans- 
verse electric and transverse magnetic modes is described resulting in an expression 
for the force on the plates in terms of the eigenvalues of the 2D Laplacian on the 
piston cross section. The force regularization is analyzed numerically and treated 
analytically resulting in absolutely convergent expressions. Asymptotic formulas are 
provided. The case of pistons with polygonal cross section is analyzed in section |4] 
and the use of the method is exemplified by computing the force on pistons with a 
triangular cross section. Finally, conclusions arc in Section [5] 

2. General formalism 

Here we summarize the formalism developed in Ref. 171 The Casimir force can be 
obtained averaging the stress tensor over the probability distribution of the quantum 
field. The stochastic quantization method of Parisi and Wul^ii^l provides a stochastic 
equation in a pseudo time, s, that reproduces the full probability distribution in 
space, r, and Wick-rotated time, r = zt, of a quantum field at temperature T. In 
the case of a scalar bosonic field of zero mass 4>, such stochastic equation reads 



where /^(t, r; s) is the noise term, mimicking the source of fluctuations for the field (j). 
It is a white noise, (5-correlated in r, r, and s, with intensity 2kBT. Equation ([T]) can 
be solved by diagonalization. Its temporal part involves the Matsubara frequencies 
uJrn — 2T:m / j3hc, m g Z, and the spatial part requires the diagonalization of the 
Laplacian in the geometry with the appropriate boundary conditions (BC), 




(1) 



Al/n(r). 



(2) 



If the stress tensor is a bilinear function of the field (as is the case in electro- 
magnetism) , T — Tr[<j>, (j), r] , which defines the stress tensor operator T, then the 
average over the fluctuations can be computed by decomposing the field (/) into its 
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eigenfunctions, to obtain (see details in Ref. |7|: 

where 7~nn(r) — T[fn, /n7^]- The sum over the Matsubara frequencies, a;^, gives, 



2 

1 



(4) 



Finally, to obtain the Casimir force over a certain body, the stress tensor must be 
integrated over the surface fl that defines the object 

Fc - </ (T(r)) • dS. (5) 
Jn 

This expression allows to compute the quantum Casimir force including the 
effects of a finite, nonvanishing temperature, in terms of the eigenvalues and eigen- 
functions of the Laplacian operator. 



2.1. Regularization 

Let us note that, as it occurs in most of the calculations of the Casimir Force, the 
expression for (T(r)) in Eq. Q is generally divergent at every point of space r if the 
sum over eigenvalues runs to infinity. Indeed, an estimation can be made considering 
Fourier modes with wavevector k. The eigenvalues of the Laplacian grow as ~ fc^ 
and the contribution of each mode to the stress tensor in the case of Dirichlet BC (as 
is the case in electromagnetism) is also Tk ~ k'^- Substitution into Eq. Q results 
in a divergent expression for the stress tensor on the object boundary. 

However, the expression for the force, that is obtained integrating over the sur- 
face of the body, Eq.([5]), is finite. It means that the integration over the body 
regularizes the divergences of the averaged stress tensor. If we assume that such 
regularization is carried out mode by mode, we can interchange the integration over 
the surface and the summation over eigenvalues, to obtain, 

i'r„„(r).dS, (6) 
Jn 

which is a finite result. Therefore, the interchange of the integral and summation 
regularizes the Casimir force, avoiding the use of ultraviolet cutoffs. Other regular- 
izations, that in some cases may lead to non-universal forces or fluctuations are, for 
instance, the subtraction of the vacuum stress tensor or by averaging the stress 
tensor over a finite area or a finite time ^^^l. Having regularized the divergences, 
Eq. ([6| , provides a new expression to calculate Casimir forces for a given geometry 
by diagonalizing the Laplace operator. So, this approach is suitable for numerical 
calculations of Casimir forces in complicated, realistic geometries. Moreover, this 
method leads to the force directly, not as a difference of the free energy with respect 
to a reference state, which in some configurations it may be difficult to establish. 



he 



'c = 



^f3hcXn 
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Fig. 1. Geometry of the considered problem. The Casimir force is evaluated for the plate, of area 
A and general geometry. The net force is obtained as the force exerted by a plate at distance L 
and another plate at distance Loo — ^ oo. The cylinder is oriented along the x axis. 

Other authors have obtained formulas for the free energy in terms of the eigenvalues 
of the spatial operator, but not for the force E^EH 

The expression for the Casimir force, Eq. (|6|, allows one to evaluate the quantum 
limit (by setting the temperature equal to zero) and the classical limit 

lim Fc = ^ V / r„„(r) • dS; lim Yc^keTy^^ I T„„(r) • dS. (7) 

The convergence of these expressions depend on the sum and integral exchange. The 
two limits, quantum and thermal, show that the driving force of the fluctuations has 
different origin. In the first case, the presence of the factor h indicates the quantum 
nature of the fluctuations, whereas in the second case, the factor ksT reveals its 
thermal origin. 



3. Electromagnetic Casimir force for a piston of arbitrary cross 
section 

3.1. Electromagnetic Casimir force 

Let us consider the geometry depicted in Fig.jl] It consists of a piston of area A, but 
general shape, made of a perfectly conducting metal surface!^. Two flat conducting 
plates of the same cross section of the piston are placed at a distance L apart along 
the X direction. The plates are perpendicular to the surface of the cylinder. We 
calculate the Casimir force, Fc, for this plate by evaluating the expression Eq. Q 
for (T(r)) and later integrating over the surface. In order to obtain a finite result, 
we need a third auxiliary plate located at infinite distance, Lao — ^ oo. 

First, we solve the eigenvalue problem for the EM field and apply Eqs. ([S]) and 
([5| to obtain the force. For the EM field the normal component of the stress tensor 
reads T^,^ = El + Bl - ^E^ - ^B^ with BC: E x n = and B • n = 0, where n 
is the surface normal vector. In this geometry, the EM field can be decomposed in 
transverse electric (TE) and transverse magnetic (TM) modes, which are discussed 
independently For the TM modes, the magnetic field is transverse to the x 
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direction, and the vector potential A can be written as 

A ^ {-CWlD,dyDda,C,d^Dd.,C)e-"^*. (8) 

Here the fields C(x) and D{r±) satisfy 

dlC{x) = -klC{x), (Neumann BC on x = 0, L) 

ViAi(r_L) = -\lDn{vA_), (Dirichlet BC on S), (9) 

where S is the surface of the cylinder and r_L = (y, z). For the TE set, the electric 
field is transverse to x, so the vector potential is A = (0, — S'SziV, S'9y-/V)e~*"*, 
where the functions S{x) and iV(rj^) satisfy Eqs. Q with the opposite BC: Dirichlet 
for S{x) and Neumann for A^(r_L). However, in this case, the constant eigcnfunction 
(with eigenvalue — 0) must be excluded as it gives A = and then E = B = 0. 

Substitution of the TE modes into the expression for the stress tensor, and inte- 
gration over one side of the plates gives, after a long but straightforward calculation, 
the Casimir force as 

1 °° 1,2 

(^""> = ^ E E E ^j^t^ta^' (10) 

where fc^ = {uxTt/L)'^. For the TM modes, one obtains exactly the same expression, 
but are the eigenvalues of the two-dimensional (2D) Laplacian with Neumann 
BC. We will denote the complete set of eigenvalues of the Laplacian with Neumann 
(excluding the zero eigenvalue) and Dirichlet BC by the index p. The expression 
above is the equivalent of Eq. ([s]) when the spectrum can be split into a longitudinal 
and transversal part, that is, A^ = fc^ + Ap. The sum over the Matsubara frequencies 
can be performed as in Q. In the limit of vanishing temperature (wi = 27r fc^T/Zic — >■ 
0) this reduces to 

hm / (nrV5. = fic7r2EE / (11) 

3.2. Regularization 



The series (11) and the one at finite temperature are divergent, but the net Casimir 
force, which is the difference between the force exerted by the plate at distance 
L, and the plate at Loo — >■ oo (see Fig. 1), is finite. To apply the regularization 
method described in Section [2T] an eigenvalue ordering method is used to obtain 
an equivalent of total force (on both sides) of each mode. That is, a kernel ordering 
method is considered. In the case of vanishing temperature this reads 



FT=n — lim hcn^ y y ni 

n^ = l p 



K 



(12) 



where K{x) is the regularizing kernel that must satisfy lim2._j.o K(x) = 1 and 
limj._5.oo K{x) = 0, and Q is the cutoff eigenvalue. 
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To study the convergence of the kernel regularization (12 1 we consider the case 



of a circular plate. In this case, the eigenvalues are the zeros of the Bessel function 
Ju{f) and its derivative. The sum can be done numerically and Fig. [2] shows Fc as 
a function of Q. It is observed that although the force on both sides diverge when 
Q — oo (see inset of Fig. [2]), the net force on the plate converges to a finite value. 




Fig. 2. Convergence of the Casimir force on a circular plate using the kernel regularization method 
with K{x) = . The circle radius is ij = 1, the interior plate distance is L = 0.5 and the exterior 
plate distance is L^o = 100.0. The n um ber of eigenvalues used in the evaluation is 3000, being 
AgQQQ ft! 77.5. Main figure: net force ( |12[ | as a function of the cutoff value Q (solid line) and the 
asymptotic value obtained using the analytic regularization \1A\ (dashed line). Inset: force on the 
interior (solid line) and exterior (dashed line) sides of the plate as a function of Q. 



For a general piston geometry, the force can be regularized analytically using 
the Chowla-Selberg summation formula to rewrite the sum over the variable in 
Eq. (10)'^. This formula extracts the divergent, L-independent part of the sum- 



mation, which cancels when the integral in Eq. ( 10 ) is performed for both sides of 



the piston. Equivalently, the L-independent contribution is the same for the second 
plate at distance L or at distance Loo, resulting in 



Fr. 



/3 



1 



(13) 



Here, A = 2tt/ (ihc is the inverse thermal wavelength. This expression gives the finite 
or regularized Casimir force between two plates at distance L, valid for any cross 
section and temperature. The precise geometry of the plates enters into the double 
set of eigenvalues of the Laplacian Ap. 



In the limit of vanishing temperature the sum over m in Eq. ( 13 ) can be rigor- 
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ously replaced by an integral, resulting in 

. oo 

Ft^o = - ^ 51 [Ko{2nLXp) + K2{2nLXp)] . (14) 

p n—l 

Here K^ix) is the modified Bessel function of order a. In Ref. TTit was obtained a 
formula for the free energy of the configuration considered here, that, after differ- 
entiation with respect to the distance between the plates, leads to the force above. 
In a similar fashion, we can calculate the thermal Casimir force when S — >■ 0. 



Then, in Eq. (13), only the term with m = is different from zero, resulting in 



3.3. Short and long distances 

For short distances, the summation over p can be calculated without explicitly 
knowing the shape of the piston and its eigenvalues. For simple cross sections, the 
eigenvalues, Ap, scale with the inverse of the square of the typical size of the piston 
because of dimensional arguments. So, for distances L much smaller than the section 
of the piston, we can replace the sum over the eigenvalues, A^, by an integral. In 
doing that step, we need the asymptotic expression for the density of states of the 
Laplacian in two dimensions, that reads (for both Dirichlet or Neumann BC^i^: 



p^{X) ^ [2XA + i^P]0{X)/47T + xS{X), 



(16) 



where A is the piston area, P its diameter, r/ — 1(— 1) for Dirichlet (Neumann) 
and X is a measure of the perimeter curvature defined as follows. For a piecewise 
continuous boundary. 



-E 

24 ^ 



TT 

ai 



TT 



(17) 



where d'jj are the line elements of the continuous segments, ^(7^) the finite curvature 
at each point of the boundary, and ai is the angle of each vertex joining contiguous 
segments. 

When integrating using the density of states, the zero eigenvalue of the Neumann 



EC must be excluded. Therefore, the expression ( 13 1 in the limit of short distances 
can be computed as 



Fc=- 



1 



poo 

E / rfAp,(A)E 

_ 1 1 "'0 



Vm2A2 + A2 



hm - E 



Vm^A^ + A2 



where the total density of states (Dirichlet plus Neumann) has been considered. 
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After a lengthy but straightforward calculation 
A 

Fc=- ksT-^ y [Li3(e-2^™^) + 2LmALi2(e-2^"^) + 2L^m^A^Lh{e-^'^"'^)] 

m—O 

oo . 

-2fcBT(2x-l)E,2zSr3I' (18) 

m— 1 

where the apostrophe in the sum indicates that the term m — O should be multiplied 
by a factor 1/2 and Lis(z) are polylog functions. Note that the force depends on 
the area and the curvature, but not on the perimeter of the piston. 
Taking the limit of low temperatures the quantum Casimir force is 

_ ficATT^ hc{2x . . 

^^=°" 240L4 24L2 • ^1^^ 

The first term is the well-known result of the EM Casimir force for infinite parallel 
plates 2] and the second term is a contribution of the curvature. 

In the thermal case, of high temperatures, the force only depends on the area, 

FZo^~^m- (20) 

In the opposite long distance limit, when L is much larger than the typical size 
of the plate only the smallest eigenvalue contributes to the sum, with the result 

Tjifar \3/2 — 2LAi T;ifar ; rrn \ — 2LAi /m\ 

2v TTi^ 

being gi the degeneracy of Ai. 



4. Polygonal pistons 

A particularly interesting case is that of pistons with a regular polygonal cross 
section. If N is the number of sides of the polygon and Rc is the exterior (circumscrit) 
radius of it, the polygon area and curvature parameter are 

sin(2^/iV) 2 N-1 

Using the value of the first eigenvalue, the far and near regimes of the force can be 
easily computed. 

To illustrate the intermediate behavior and the transition from the near to the far 
regime, we study the case of a equilateral triangle of side a. The relevant quantities 
for the computation in the asymptotic regimes are'^ 

A = a^VS/A; x = 1/3; >h = IQt^^ /^a^ (Neumann); gi = 2. (23) 

At intermediate distances, that is, L comparable to the size of the plates, one must 
solve the eigenvalue problem which can be obtained numerically using finite element 
methods for a general geometry or analytically for special geometries as the present 
casein. The dependence with the side of the triangle is \p{a) — Ap(a = l)/a^, so in 
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the quantum limit the Casimir force as expressed in Eq. (14) is a function of L/a 
when the force is multiplied by , while in the thermal limit the Casimir force Eq. 
( 15 1 is a function of L/a when the force is multiplied by a. 

Figure |3] (left panel), shows the Casimir force in the quantum limit, Eq. ( [l4| , 
for an equilateral triangle of side a. To compare, five results are presented: (1) the 
full computation using an extensive list of eigenvalues {N = 4000, that ensures con- 
vergence for L/a > 0.03) that we consider as the true value; (2, 3) the asymptotic 
expression for near plates ([19^ with and without the curvature correction; (4) the 



asymptotic expression for far plates ( 21 ); (5) a partial sum considering only the first 
20 eigenvalues. A similar comparison in the thermal limit is presented in the same 
figure (right panel). The analysis of both cases, quantum and thermal, is qualita- 
tively similar. When the plates are separated, only a small number of eigenvalues 
is needed to obtain the force (20 eigenvalues for L/a > 0.2). At smaller distances, 
more eigenvalues would be needed, but there is an excellent match with the asymp- 
totic expansion at short distances. The curvature correction in the quantum limit 
allows to extend its validity up to L/a « 0.4. 




Fig. 3. Left: Quantum limit of the Casimir force for a piston of triangular cross section of side a 
as a function of the distance L between the plates. Solid circles are obtained summing over 4000 
eigenvalues of the Laplacian, using Eq. |l4]l. Open squares (overlapping with circles for L/a > 0.2) 
are obtained summing only over 20 eigenvalues. Asymptotic expressions for near plates \19\ without 
the curvature correction (dotted line), with the curvature correction (solid line) and the asymptotic 
expression for far plates | |21| (dashed line). Right: Thermal limit of the Casimir force for the same 
configuration. Solid circles are obtained summing over 4000 eigenvalues of the Laplacian, using 
Eq. l |15| . Open squares are obtained summing only over 40 eigenvalues. Asymptotic expressions 
for near plates \20\ (solid line) and the asymptotic expression for far plates \21\ (dashed line). 



5. Conclusions 

The new method described in Ref. [7|to compute Casimir forces using the stochastic 
quantization formalism has been used to calculate the force on pistons of arbitrary 
cross section. The method is quite simple and it provides the force directly, instead 
of the free energy. It only requires spectral decomposition of the Laplacian opera- 
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tor in the given geometry, summation of the eigenvalues and the integration of the 
eigenfunctions along the boundary of the object. Such integration of the eigenfunc- 
tions over the surface of the body leads to a regularization of the Casimir force, 
producing finite results for the force. Quantum and classical limits are recovered as 
well as intermediate results for finite temperatures. The electromagnetic Casimir 
force on pistons with arbitrary cross section is analyzed in detail. The computa- 
tion of the 2D eigenvalue problem allows to obtain the Casimir force in the whole 
range of separation distances. To demonstrate the method, the force on a piston of 
triangular cross section is computed numerically. 
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